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We generalize the pointwise decay estimates for large data solutions of the defocusing 
semilinear wave equations which we obtained earlier under restriction to spherical symmetry. 
Without the symmetry the conformal transformation we use provides only a weak decay. This 
can, however, in the next step be improved to the optimal decay estimate suggested by the 
radial case and small data results. This is the first result of that kind. 



I. INTRODUCTION 

This Article builds up on our previous result [BS09] where we considered wave equations with 
the defocusing nonlinearity 

□</.=-i0rv, = 92 -A, (1) 

in n = 3 spatial dimensions restricted to spherical symmetry. Here, by generalizing the technique 
we can remove the symmetry condition and prove the same decay result 

C 

\4^{t^x)\<- ^ (2) 

' - {l + t+\x\){l + t-\x\)P-^ ^' 

with some constant C depending only on initial data and the power 3 < p < 5. For the initial data 
((/>0,</>i) we assume the regularity 

(0,5t0)t=i = {<t)Q,4>i) e C3(M3) X (3) 

implying existence of classical solutions cf) G C2(M x M^), and compact support"*^ in \x\ < a := ^ 
which, among others, guarantees that the positive definite energy 

E = I {\\dt4>? + \\y<t^? + -^Mr'] d'x (4) 

is finite but not necessarily small. 

Under these conditions Jorgens [Jor61] has shown global existence, Pecher [Pec74] uniform 
boundedness and Strauss [Str68] weak uniform decay More references to similar results can 

be found in our previous publication on that equation [BS09]. We are not aware of any stronger 
pointwise estimates existing in the literature. While under spherical symmetry the decay rates 
were well motivated by the numerical analysis [BC08] and an asymptotic result for small initial 
data [SBCR07] stating that for large t and fixed r the solution behaves like 

,^(t,r) = ^ + 0(t-^), 



^ Symmetries of the equation (1) allow for mapping any compactly supported initial data onto the region \x\ G [0, a[ 
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such strong statements are lacking in the non-spherical case. We are only guided by our spherical 
result for big data [BS09] and a series of results giving pointwise decay estimates for small initial 
data [Asa86, ST97, Szp08]. In this context the estimate (2) seems optimal. 

Here, we combine the technique of conformal compactification developed by Choquet-Bruhat, 
Christodoulou and others in [CBPS83, Chr86, BSZ90], the boundedness result for big data of 
Pechcr [Pec 74] and some pointwise estimates developed by Asakura [Asa86] and one of us [Szp07]. 

Our method has a few steps. First we perform a conformal compactification of the space-time. 
We map the double-null coordinates u = t + r, v = t — r where r = |x| and the angular coordinates 
6, ip according to 

u := , V := — , 6 -.= 6, S := ip. (5) 

u V 

The conformal factor := r/r multiplies the transformed solution (f>(u,v) := ^~^(f){u,v) which 
satisfies the transformed wave equation 

+ {uv)P-^^\^\P-^ = (6) 

in a precompact region of spacetime. There, we are able to show boundedness of some pseudo- 
energy flux which we further use to show the uniform boundedness of |(/>(S, u)| < C. Finally, the 
inverse transformation gives us the weak pointwise estimate 

\(f)iu,v)\<C -^1 = ^. (7) 

In the radial case (cf. [BS09]) there was a family of conformal transformations available, being the 
mappings of the pair {u,v), which could be adjusted to the power of the nonlinearity p in such a 
way that the conformal factor 17 already delivered the desired optimal pointwise decay. Without 
the spherical symmetry these transformations are no more conformal except the one which was 
related to p = 3 and is now being used. Hence, the rate of decay in (7) corresponds to that of the 
radial problem with p = 3. 

In order to improve it (for p > 3 the decay is faster), we need a second step in which we make 
use of an integral representation of (1), known as the Duhamel formula, and its estimates, briefly 

II. CONFORMAL TRANSFORMATION 

On R X M3, let r+ := {{t,x) | < |x| < t} and T" := | < |x| < -t} denote the 

interior of the forward and backward lightcone of the origin respectively, where the coordinates 
i : M X — )• M and a; : M x — )• M"^ are the canonical projections. Consider the map 

$:r+^r-, {t,x)^th^, (9) 

— 

It is an analytic bijection with analytic inverse 

{-t,x) 



f^-x 
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If 77 := di^ — Sij dx^dx^ denotes the Minkowski metric on R x R^, then its pullback ^*r] by $ 
satisfies 
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on T"*", where 



as ^*dt = {2t^n - l)n dt - 2txin'^dx' and ^*dx' = -2tx'n^ dt + {5] + 2x'xj^)ndxK Hence, 
the mapping # is also conformal. The conformal factor Q, is analytic and has the property that 
= Q.\x\. 

Prom the calculations in the spherically symmetric case, more precisely from Equation (10) in 
[BS09], it easily follows that if /i € C'^{T~) is a twice continuously differentiable function on T~, 
its pullback G C'^{T^) and satisfies 

U{rL<^*h) = Q!^<^*Uh. (10) 

In particular, let G C'^iT^) be a classical solution of the semilinear wave equation 

□(^ + (^101^-^ = (11) 

for some p > 2, then its conformal transformation := $*(f2~^(/>) G C'^{T~) satisfies 

+ ($*j^)p-vi^r^ = 0, (12) 

where denotes the push-forward by the diffeomorphism On the other hand, it follows directly 
from (10) that if -0 G C'^{T~) is a classical solution of the transformed equation (12) on T~, the 
function (f) := Q^*ip G C'^(T^) solves the original equation (11) on T^- Note that for p > 3, the 
non-linearity in Equation (12) is regular at the origin. In the case p = 3 the functions cj) and i/' 
solve the same Equation (11), which is in this case called conformally invariant. 

The method of improving a boundedness result for into a decay estimate for ^ by means 
of the conformal transformation described above is directly related to the Morawetz vector field 
[Mor62] 

Z := u^du + v'^dy (13) 

with null coordinates u := t + \x\ and v := t — \x\ used in the context of vector field methods to 
obtain decay. In fact, Z is the pullback by $ of the time-like Killing vector field dt on T~, 

Z = {t'^ + x'^)dt + 2tx-V = ^*dt. 



III. LOCAL ENERGY ESTIMATE 



Having in mind the transformed wave eqiiation (12), the aim of this section is to establish a 
local energy estimate for equations of the form 

D'ilj + c{t,x)i;\'ilj\P-^ =0 (14) 

with a non-negative function c of class C^. An energy density which is expected to be useful for 
such equations is 

e:=l{dt^r+^mf+^c\^i;r\ 

While the associated energy will, in general, not be conserved, it will be sufficient to prove bound- 
edness of on the relevant region of T~ , provided the function c is monotonically decreasing in t 
there. 
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Let V' € C^(T ) be a classical solution of (14) on T for p > 2. Then the vector field S given 



by 



is continuously differentiable and 



div£^ 



S := edt — (ptil^ V'^) • grad 



p + 1 



{dtcM 



P+i 



p + l 



holds on T~ ■ Assume furthermore that dtc < is non-positive then the same is true for divS. 
This implies, recalling the assumption that the initial data is compactly supported in an open ball 
of radius a = 1/2 about the origin, that 



•Di(-1,0) 



> / (div£")a; = d{ 
Js Js 

= / {i£dt)dx^ A dx^ A dx^ — / i£LU 

iDi(-l,0) J^Dc{l,0) 



*i3a(l,0) 



L>i(-1,0) 



(1,0) 

P+i 



D„{1,0) 



where denotes the interior multiplication with £, uj = dtAdx^ Adx'^ Adx^ is the standard volume 
form and all hypersurfaces shall be oriented by the transverse vector field dt- The region S C T~ is 
defined as this part of the future of the image $£)«(!, 0) which lies in the past of the hypersurface 
{t = —1}. Here, Dp{t,x) := {{t,y) \ \y — x\ < p} is an open disk of radius p around {t,x). As the 
integral over Da{l, 0) depends only on the initial data smd (pi there is a constant C such that 



Oa(l,0) 



< C. 



Therefore it follows that 

7Di(-1,0) 



1 



1 



< c. 



(15) 



The energy Eq now controls past light-cone integrals of cl^l^"*"^ or, more generally, of the flux 
density 

2 



y — X „ 
dt + -f T ■ V 

\y-x\ 



If y — X 



\y - x\ 



1 



p + l 



c|V| 



at any point {t,x) £ Q := H { — 1 < t < 0} according to the following Proposition 1. For 
notational convenience, let K{t,x) := {(s,y) | — 1 < s < t, \y — x\ < t — s} he the solid truncated 
backward light-cone at (t, x) and M{t,x) := {(s,y) | — 1 < s < t, \y — x\ = t — s} \ts mantle, 
where {t,x) is a point in the relevant region Q. 

Proposition 1. Let G C'^{Q) he a solution of (14) with initial data supported in Di(— 1,0) 
satisfying the estimate (15) . // dtc < on Q then 



Flux(t, x) :-- 



V2jM{t,x)['^\\y-x 



If y — X 



dttp 



+ 



1 



P + 



P+l 



< En 



holds for any {t, x) E Q. 
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Proof. Fix {t,x) G Q. Since S is continuously difFerentiable, direct application of Stokes's theorem 
yields 



/ divS = / igu! — f 

JK{t,x) JM(t,x) JDi+ti-i 

1 



V2 



M(t,x) 



{-l,x) 



\y - x\ 



I 

JDi+t{-l,x) 



isdt 



If y — X 



2\\y-x\ 



\/2 JM{t,x) 

= Flux(i, x)- 

JDi+t{-l,x) 

Hence, it follows that 

Flux(i,a;) = / e + / divS 

JDi+ti-l^x) Jk(Ux) 



Di+t(-l,a;) 



<-L 



'Di(-1,0) 



+ 



divf 



K{t,x) 



□ 



IV. BOUNDEDNESS AND WEAK DECAY 



Using the local energy estimate established in Proposition 1 the boundedness of ^ in the region 
Q = 7 n{— l<t<0}is quite immediate. The argument is exactly the same as in [BS09] and 
goes back to Pecher [Pec74]. 

Theorem 1. Let ip € C'^{Q) he a solution of (14) with 2 < p < 5 and initial data supported in 
1,0) satisfying the estimate (15). If c > is uniformly bounded and dtc < on Q then ^ is 
uniformly bounded. 

Proof. Let ipo G {Q) be a classical solution of the homogeneous equation Oipo = with the same 
initial data as ip. Then, for a fixed {t,x) G Q, it holds that 



c{s,yM\^is,y) 



dBt-s{x) 



t - S 



dy ds. 



Due to the fact that 2 < p < 5 there exists a q with 3/2 < q < {p + l)/{p— 1). Changing variables 
and applying Holder's inequality yields 



/'/ 

J -I Jd 



dy ds 



dBt-six) 



t-s 

c\ip\P{t - \y-x\,y) 



Bi+t(x) 



\y - x\ 



dy 



(16) 



< / c'?|Vr(t-|y-x|,y)dy 

UBi+t(x) 



lJBi+t{x) 



\y-x\ 



dy 
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Consider the first integral. Since q> 1 and c is bounded, so is ^. Furthermore, < pq— (p+1) < 
q, so that 



ci\i;r{t-\y-x\,y)dy 



Bi+t{x) 

<c( sup ||V.(r,.)||L^y'^ / clV^rS 
where < 7 < 1 is such that 79 = M — (p + !)• By virtue of Proposition 1 

^ / cW^'^ <{p+l) Flux(t, x) < CEo. 

V2 JM(t,x) 



The second integral in (16) can be calculated directly to give 



y-a;| "-^ dy = An r "-^ dr = Air -{I + t)^"-^ < C 

Jo 



2q-3 



'Bi+t{x) 

because q > 3/2. To sum up, the estimate 



\i^-i^o\{t,x)<CE^i sup mT,-)\\Lo 

\-l<T<t 

holds true for any (t, x) G Q. But since the solution of the homogeneous equation is clearly 
bounded and < 7 < 1 this estimate implies the boundedness of itself uniformly on Q. □ 

As a short digression, note that the proof of Theorem 1 given above uses the boundedness of 
light-cone integrals of the quantity cl-^j^^^ which is controlled by the "potential" part of the flux 
defined in Proposition 1. Alternatively, following Shatah and Struwe [SS98], the integral in (16) 
could have been estimated by 



c\ip\^it -\y- x\,y) 



Bi+tix) 

< 



\y - x\ 



dy 



1 

/ c^\^pf(P-'\t-\y\,x + y)dy " [ 

JBi+tiO) J Ue 



\'ip\^{t-\y\,x + y) 
/Bi+t(o) |yP 



dy 



(17) 



The first integral would then be controlled by 

i^+trm&[Mlx)) 

when 1 < p < 4 and by 

when 4 < p < 5. But now, contrarily, the "kinetic" part of the flux could be used to estimate 
both the L^-norm of on M{t, x) by Sobolev embedding as well as the second integral in (17) by 
Hardy's inequality. 

With the function ^ bounded on Q a decay of (p towards the future follows directly. 
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Corollary 1. Let (/) be a classical solution of the wave equation (11) for 3 < p < 5 with initial data 
(j)Q G C^(M'^) and (pi G C^(M^) given at t = 1 and which exists globally towards the future. Assume 
that the support of cpo o-nd 4>i is contained within the open ball Ba{0) of radius a about the origin. 
Then there is a constant C > such that (j) satisfies the decay estimate 



l'^(^'"^l^ (i + t + N)(i + t-|x 



^ (18) 



for all t > 1 and x G 



Proof. Given such a solution cp of class C^, it was shown in Section II that its conformal transfor- 
mation ip = $*(ri~^0) is a classical solution of the wave equation (14) on the future of $L'q,(1,0) 
in T~ with 

C = = (t2 _ a;2^P-3 

according to equation (12). So c is certainly bounded on the future of $Da(l,0) in T~ because 
p > 3 and 

dt [^*np-^] = 2{p - 3)^$*^^f-^ 

which implies dtc < on again by reason of p > 3. Moreover, as detailed in Section III, the 
assumptions on the support of cpQ and (pi guarantee that the support of the transformed solution 
restricted to {t = —1} is compactly contained in Di(— 1,0) and that therefore the estimate (15) 
holds. Thus, Theorem 1 applies and yields boundedness of ip on Q. Since ip is also bounded on 
the compact region S it follows that ip is bounded on the whole future of $Dq.(1,0) in T~, say 
IV'I < C. But then 

\(t){t,x)\ = n{t,x)\^*tp{t,x)\ < cn{t,x) 

on the whole future of Da{l,0) in T"*". In this region, t + |x| > 1 and i — |a;| > 1 — a, so that there 

1 2(1 + ^ 

^(*'^) - ^23^ ^ ^i + t + \x\){i + t-\x\y 

Hence, the bound of \<p\ can be written in the regularized form (18). □ 

V. IMPROVEMENT OF THE DECAY ESTIMATE 

The pointwise decay result of Corollary 1, 

can be improved by applying this bound to the nonlinear term in the wave equation (1) and solving 
the wave equation by inverting the wave operator, i.e. (j) = {—\(p\P~^ cp) + X(t>o,(t>i- Here, X^c^i 
represents the contribution from the initial data (3), i.e. solves the linear wave equation = 
with {Xjdtx)\t=i = (00) 0i) ^-iid is well-known to decay as |x(t,a:)| < C/t. Due to the Huygens 
principle in three dimensions xit, x) is supported in the outgoing light-cone 1 — a < t — \x\ < 1 + a. 
Hence, its decay can be written as 

\xit,x)\<7- r-FT7 r-F^ (20) 

''^^ ' - {l + t+\x\){l + t-\x\)i ^ ' 
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with any power q and some C depending on q. Since ^ G is a classical solution and □ ^ is a 
positive measure on M x we can estimate 

wt, x)\ < + \x(t, x)\ < a-^- ---^ + --^ (21) 

'^^ ' ^' - '^^ ' ^' - {l + t + \x\)P{l + t-\x\)P {l + t + \x\){l + t-\x\)i ^ ' 

The inverse wave operator can be represented by the Duhamel integral formula and bounded 
pointwise. According to the Lemma 1 from [Szp07] we get for p> 2 

□-1 \ < ^ (22\ 

{l + t + \x\)P{l + t - \x\)v - (1 + i + \x\){l + t - \x\)P-^ ■ ^ ^ 

Choosing q = p — 2 we arrive at our main result 

Corollary 2. Under the assumptions of Corollary 1 there is a constant C > such that (j) satisfies 
the improved decay estimate 

for all t >1 and x G . 



VI. OUTLOOK 



In some sense, a very similar problem, the linear wave equation with a strong positive potential 

d^^ -A(j) + V{x)(l) = 

having prescribed decay at spatial infinity V{x) ^ l/|a;|^ still lacks a sharp pointwise decay esti- 
mate. It also has a positive definite energy and can be conformally transformed to a form analogous 
to (14). However, the function c{t,x) is no more regular at t = \x\ = 0. Nevertheless, we expect 
that our method can be extended to cover this weakly singular case, too. 



Acknowledgments 

The authors want to express their gratitude to the Mathematical Institute in Oberwolfach and 
to the organizers of the workshop "Mathematical Aspects of General Relativity" (Fall 2009) for 
hospitality and creative atmosphere during the work on that project. One of the authors (NS) 
also wants to thank the organizers of the workshop "Quantitative Studies of Nonlinear Wave 
Phenomena" at the Erwin Schrodinger International Institute for Mathematical Physics in Vienna 
(Winter 2010) during which this work has been finished. 



[Asa86] F. Asakura. Existence of a global solution to a semi- linear wave equation with slowly decreasing 

initial data in three space dimcnstions. Comm. Part. Diff. Eq., 13(11):1459-1487, 1986. 
[BC08] P. Bizon and T. Chmaj. 2008. private communications. 

[BS09] R. Bieli and N. Szpak. Global pointwise decay estimates for defocusing radial nonlinear wave equa- 
tions. 2009. arXiv: 0903.0799 [math.AP]. 

[BSZ90] J. C. Baez, I. E. Segal, and Z.-F. Zhou. The global Goursat problem and scattering for nonlinear 
wave equations. J. Func. Anal, 93:239-269, 1990. 



9 



[CBPS83] Y. Choquet-Bruhat, S. M. Paneitz, and I. E. Segal. The Yang-Mills equations on the universal 

cosmos. J. Func. Anal, 53:112-150, 1983. 
[Chr86] D. Christodoulou. Global solutions of nonlinear hyperbolic equations for small initial data. Comm. 

Pure Appl. Math., 39:267 282, 1986. 
[Jor61] K. Jorgens. Das Anfangswertproblem im Grossen fiir eine Klasse nichtlinearer Wellengleichungen. 

Math. Z., 77:295-307, 1961. 
[Mor62] C. S. Morawctz. The limiting amplitude principle. Comm. Pure Appl. Math., 15:349-361, 1962. 
[Pec74] H. Pecher. Das verhalten globaler Losungen nichtlinearer Wellengleichungen fiir grofie Zeiten. Math. 

Z., 136:67-92, 1974. 

[SBCR07] N. Szpak, P. Bizon, T. Chmaj, and A. Rostworowski. Linear and nonlinear tails II: spherical sym- 
metry. 2007. arXiv: math-ph/0712. 0493; accepted for publication in Journal of Hyperbolic Differential 
Equations (JHDE). 

[SS98] J. Shatah and M. Struwc. Geometric wave equations, volume 2 of Courant Lecture Notes in Mathe- 
matics. New York University Courant Institute of Mathematical Sciences, New York, 1998. 

[ST97] W. Strauss and K. Tsutaya. Existence and blow up of small amplitude nonlinear waves with a 
negative potential. Discr. Cont. Dynamical Systems, 3(2):175-188, 1997. 

[Str68] W. Strauss. Decay and asymptotics for Du = F{u). J. Func. Anal, 2(4):409-457, 1968. 

[Szp07] N. Szpak. Simple proof of a useful pointwise estimate for the wave equation. 2007. arXiv: math- 
ph/0708.2801. 

[Szp08] N. Szpak. Linear and nonlinear tails I: general results and perturbation theory. Journal of Hyperbolic 
Differential Equations (JHDE), 5(4):741-765, 2008. arXiv: math-ph/0710.1782. 



